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Abstract 

We present the results for three-loop beta-function for the Higgs self-coupling 
calculated within the unbroken phase of the Standard Model. We also pro- 
vide the expression for three-loop beta-function of the Higgs mass parameter, 
which is obtained as a by-product of our main calculation. Our results coin- 
cide with that of recent paper arXiv:1303.2890, In addition, the expression 
for the Higgs field anomalous dimension is given. 



The Higgs self-interaction coupling being the fundamental parameters of 
the Standard Model (SM) Lagrangian describes the interactions of Higgs field 
with itself and is strongly related to the Higgs mass via electroweak symmetry 
breaking. Having in mind the discovery of the Higgs boson [H [2] the Higgs 
self-interaction coupling can be deduced directly from the experimental data. 
In order to obtain a very precise SM prediction for the running Higgs self- 
coupling at some high energy scale, one usually uses value extracted from 
Higgs mass measurements around electroweak Mz scale. Objects of interest 
are scales up to the Planck mass, so one inevitably makes use of renormaliza- 
tion group equations (RGE) to connect these scales. The SM parameters in 
such a studies are usually defined in the minimal subtraction (MS) scheme, 
in which counter-terms and beta-functions have a very simple polynomial 
structure. One can use RGE for finding the scale where New Physics should 
enter the game, e.g., to unify the interactions or stabilize the Higgs poten- 
tial laaiEiEiEiiH]. 

One- and two-loop results for SM beta functions have been known for 
quite a long time P[ini[IIl[Il[I3l[ia[I3[IS[ni[ISl[ISl[2ni[2Il[22]andare 
summarized in [23] . 

Not long ago full three-loop beta-functions for gauge couplings [211 [25] and 
Yukawa couplings [26] were calculated. The beta-functions for the Higgs self- 
coupling and top Yukawa coupling were also considered at three loops [27] . 
However, in Ref. [2^ all the electroweak couplings were neglected together 
with the Yukawa couplings of other SM fermionsjj 

In this paper, we provide the full analytical result for the three-loop beta- 
functions of the Higgs-self coupling A and the Higgs mass parameter m^. We 
take into account all the interactions of the SM, restricting Yukawa sector to 
include only the heaviest fermion generation. 

Let us briefly recall our notation. The full Lagrangian of the unbroken SM 
which was used in this calculation is given in our previous paper [25J. How- 
ever, we do not keep the full flavor structure of Yukawa interactions but use 
the following simple Lagrangian which describes fermion-Higgs interactions 
and the Higgs field self-interaction 

-^Yukawa = -yt{Q'^^)tR - yb{.Q^)hR - yr{L^)TR + h.C. , (1) 

£h = {D^^)^ - , (2) 

VHm = \{^^^f = x{^-^ + <P^r'^ , (3) 

^During the preparation of this paper, the authors of Ref. extended their resuh and 
incorporate the dependence |28] on the electroweak gauge couphngs and Yukawa couphngs. 
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with Q = {t,b)L, and L = {i'tit)l being SU(2) doublets of left-handed 
fermions of the third generation, ur, tpt, and r^j are the corresponding right- 
handed counter parts. The Higgs doublet $ with Yw = 1 has the following 
decomposition in terms of the component fields: 

Here a charge-conjugated Higgs doublet is introduced ^'^ with Yw = — 1- The 
Higgs self-coupling A entering tree-level Higgs potential ([3]) is of our primary 
interest. We do not add a quadratic (mass) term m^$"''$ to the potential Vh, 
since the running of the mass parameter can be deduced by considering 
renormalization of composite operator The treatment is essentially the 
same as in Ref. [27]. Some details can be found below. 

For loop calculations it is convenient to define the following quantities: 

gl gl gl v1 vl vl ^ ^ ^ ^ A 



3 167r2 ' 167r2 ' levr^ ' IOtt^ ' lOvr^ ' 167r2 ' 167r2 

where we use the SU(5) normalization of the U(l) gauge coupling g^. We 
also stress that the calculation is carried out in a general linear gauge, in 
which the vector boson propagators has the form 
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^Q = l-^Q- (6) 



A minimal way to test gauge invariance at the end of calculation is to keep 
at most a single power of ^g, which corresponds to a first order expansion of 
the result around the Feynman gauge. 

The A beta-function is extracted from the corresponding renormalization 
constant which relates the bare coupling to the renormalized one in the MS- 
scheme. The latter can be found, for example, with the help of the following 
formulae: 

7,,., 7 Z, , ,^ ,_ 

(7) 



Zhhhh _ ^xxxx _ Zhh4,+4,- 



Z\ ^/ia/ Z^+Z^~ ' 

where Zhhhh-, -^xxxx' Zhh<i>+<i>- are the renormalization constants for the four- 
point vertices involving four components of the Higgs doublet 

The renormalization constant Zh = Z^ = Z^± = Z$ can be found from 
the corresponding self-energy diagrams. It turns out that due to the gauge 
symmetry all Higgs doublet components renormalize in the same way. More- 
over, the same reasoning can be applied to the considered Higgs vertices 
giving, e.g., Zhhhh = ^xxxx = Zhh<f,+^-- 
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h 0- h .>^(f)~ \(f)- h .>^X \(l>- h .>^h \ (f)- 

a) b) c) d) 

Figure 1: For the calculation of the beta-function of the Higgs self-coupling 
we evaluate all diagrams included in Fig. a. For the calculation of the beta- 
function of the Higgs mass parameter we should insert the mass operator 
into self-energy diagrams for Higgs fields. Effectively the mass operator is 
equivalent to the quartic Higgs vertices with two external Higgs fields 
All diagrams corresponding to Fig. h should be multiplied by the factor 1/2. 



Renormalization constant for the Higgs mass parameter can be easily 
extracted from the calculations of the renormalization of the Higgs self- 
coupling. In the most simple way this can be done from our calculation 
of hh(j)'^(f)~ vertex. For this purpose we have labeled all quartic Higgs vertices 
and have extracted the results, which contain such vertices with external 
fields. This trick is illustrated on Fig. [TJ Effectively, this trick is equiv- 
alent to the insertion of the local operator = as in Ref. [27]. The 
corresponding renormalization constant can be extracted from our results in 
the following way: 



z 



hh [(/>+(/>- 



In order to extract a three-loop contribution to the considered renormaliza- 
tion constants, it is sufficient to know the two-loop results for the gauge and 
Yukawa couplings and the two-loop expression for the Higgs self-interaction. 

The relation between the bare and renormalized parameters can be writ- 
ten in the following way 

afc.Bare/i"^^'" = Za^aM = Ofc + 4"^— , (9) 

n=l ^ 

where pk = I for the gauge and Yukawa constants, pk = 2 for the scalar 
quartic coupling constant, and pk = for the gauge fixing parameters. The 
bare couplings are defined within the dimensionally regularized [22] theory 
with D = 4 — 2e. The four-dimensional beta-functions, denoted by (3i, are 
defined via 

dai{p,e) 



dhi p^ 



A = + P?^ + /3f + • • • (10) 
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1-loop 


2-loop 


3-loop 


hhhh 


246 


40 905 


8 659436 


xxxx 


246 


40 905 


8 659436 


hh(f)+(f)- 


146 


29 289 


6 741 584 




146 


29 289 


6 741 584 


0+0-0+0- 


193 


35 211 


7 597252 


hhxx 


168 


32 469 


7378 694 


hh 


9 


419 


41 369 


XX 


9 


419 


41369 


0+0- 


8 


394 


39122 


hh [0+0-] 


9 


900 


140 979 



Table 1: The number of diagrams for calculations up to three- loop order. 



with Pi being the /-loop contribution to the beta-function for aj. The ex- 
pression for Pi can be extracted from the corresponding renormahzation con- 
stants with the help of 

Here, again, stands for both the gauge couplings and the gauge-fixing. 

As in our previous work [301 ED 122] all calculations were performed with 
FORM [33], using FORM package COLOR [S] for evaluation of the color 
traces. Feynman integrals are evaluated by the method from Refs. [351 [36] 
and our own implementation of the Laporta's algorithm [37] in the form 
of the MATHEMATICA package BAMBA with the master integrals from 
Ref. [M]. 

According to the prescription of Refs. [351 ES], we introduce an auxiliary 
mass parameter M in all propagator denominators and perform an expansion 
in external momenta. Due to this, we only have to consider vacuum integrals 
with one mass scale. The subtlety of the method is related to the fact that 
one needs to introduce mass counter-terms for all the boson fields of the 
model, i.e., gauge and Higgs fields. Moreover, we have to consider diagrams 
with counter-term insertions for all the vertices of the model. 

For dealing with a huge number of diagrams (see Table [1]) instead of Fey- 
nArts [39j package, exploited in our recent studies [251 [26], we use a program 
DIANA ^0], which calls QGRAF [H] to generate all diagrams. By means 
of a prepared script we map topologies generated by DIANA on previously 
defined auxiliary topology for IBP identities. The model file for the unbro- 
ken SM, used previously with FeynArts, was converted to the DIANA model 
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format, which allow us to introduce vertex counter-terms in a convenient 
way. Most of these counter-terms were generated from the known two-loop 
renormalization constants for SM parameters and fields. We have calculated 
XXXX cLnd hh(f)'^(j)~ vertices and the obtained results are the same. The last 
vertex was used for the calculation of the /3-function for the mass parameter. 
Both calculations were performed in the linear gauge (see Eq. (^) and we 
keep only the first power of ^q. These results can be found online as ancillary 
files of the arXiv version of the paper. Since we are considering three-loop 
Green function with scalar external legs at zero external momenta it is easy 
to convince oneself (see, e.g., reasoning given in Refs. |l2l EH [251 EH]) that 
the naive anticommuting prescription for the 75 matrix is sufficient for our 
current study. 

As a result of our calculation we obtain the expressions for the three-loop 
Higgs self-coupling (|T^ and mass parameter (ITTI) beta-functions (A = ax): 



n{i) To\2 9aiA 9a2\ 27af 9aia2 9a^ 
Px ~ -'-2A 77; ;:; 1 — H — 1 



2 
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with uq corresponding to the number of SM generations and ^3 = C(3)- To 
save space we substitute all the color invariants by the corresponding values 
{Cp = 4/3, Nc = 3,Ca = 3). These results with all color invariants can be 
found online as ancillary files of the arXiv version of the paper. 

In addition, we present the expression for the Higgs field anomalous di- 
mension in the Landau gauge, which is usually adopted in the effective po- 
tential calculations HSl Hiff 
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^The result in a general linear gauge can be found online as ancillary file of the 
arXiv version of the paper. 
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27 27 27 3 831ogot 8310^0^ 

— —O'lO'TO't —0'20'TO't —0'20'bO'T — -O-bO-rO-t H 77 1 77 

10 2 2 2 lb Id 

135otA2 135o;,A2 45o^A2 .r 2t .r x2 
^ — + 45ot^A + 45o^A + 4502A2 

+ 180^04 + 180^0(2 + 18ogo^ + 18ofeo2 + ISo^A + 9oiP - 36A^ . (20) 



It should be noted that the expressions for Px and are free from 
gauge-fixing parameters ^g,^w and C,b which are present in the renormahza- 
tion constants for the considered Green functions. The one- and two-loop 
corrections are in a full agreement with Refs. [TSl ESI 112] • The contribu- 
tions f lT^ and f[T7|l coincide with the result of Refs. [2S]- In the limit of 
vanishing coupling constants 01,02,06, and Qr, the result for the Higgs field 
anomalous dimension coincides with the expression presented in Ref. [2Zj • 

To conclude, in this paper we present the expressions for three-loop renor- 
malization group quantities of the SM Higgs sector, i.e., f3\, f3m2. Moreover, 
we provide the result for the anomalous dimension for the Higgs field. The 
former can be used in a study of high energy behaivoir of the SM parameters. 
The latter may be exploited in a more accurate analysis of the Higgs effective 
potential (see, e.g., Ref. [I5lll6]). 
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